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SUMMARY 


A  calculation  method  for  transonic  separated  flow  about  two-dimensional 
airfoils  at  incidence  is  presented  in  this  report.  The  method  is  capable 
of  predicting  the  effects  of  both  leading-  and  trai ling-edge  separated  flows, 
although  the  former  is  usually  associated  with  a  complete  collapse  of  the  air¬ 
foil  flow  field,  in  which  case  the  flow  is  no  longer  supercritical.  A  viscous 
potential  flow  iteration  procedure  provides  the  connection  between  potential 
flow,  boundary  layer  and  wake  modules.  The  separated  wake  is  modelled  in  the 
potential  flow  analysis  by  thin  sheets  across  which  exists  a  jump  in  velocity 
potential.  These  sheets  are  analogous  to  vorticity  sheets  in  incompressible 
flow.  The  basic  potential  flow  method  is  a  modification  of  Jameson's  full 
potential  method.  Calculations  for  four  different  airfoils  have  been  compared 
with  experiment  for  pressure  distributions  as  well  as  integrated  forces.  The 
agreement  between  theory  and  experiment  is  generally  good  even  when  shock 
waves  are  present. 
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1.0  INTRODUCTION 


A  calculation  method  for  transonic  separated  flow  about  two-dimensional 
airfoils  at  incidence  is  presented  in  this  report.  This  aspect  of  aero¬ 
dynamics  has  had  relatively  little  attention,  although  its  importance  cannot 
be  overstated.  Separated  flow  is  present  in  almost  all  applications  of 
Dractical  interest,  and  a  full  analysis  has  not  yet  been  very  successful. 

This  is  indeed  a  challenging  problem  that  requires  not  only  a  good  inviscid 
transonic  flow  calculation  method,  but  also  an  accurate  prediction  of  boundary 
layer  development.  Above  all,  in  order  to  obtain  efficiently  any  realistic 
solutions  for  which  the  large  separated  region  exists,  the  wake  region  must 
be  modelled  appropriately. 

Only  a  few  investigators  have  attempted  to  model  this  complicated  flow. 
Hicks  (Ref.  1)  has  coupled  an  optimization  technique  to  an  existing  two- 
dimensional  transonic  airfoil  code.  He  has  been  able  to  obtain  separation 
profiles  where  the  pressure  remains  constant  along  the  separating  streamline. 
Unfortunately,  this  procedure  has  not  yet  resulted  in  good  agreement  with 
experimental  data.  Diewert  (Ref.  2)  has  treated  the  two-dimensional  airfoil 
with  separation  using  a  time-dependent  Navier-Stokes  code.  Initial  results 
for  an  18%  thick  circular  arc  airfoil  with  a  very  small  amount  of  trai ling- 
edge  separation  are  quite  qood;  however,  as  the  separated  region  increases, 
agreement  with  experiment  deteriorates.  Recently,  Barnwell  (Ref.  3)  proposed 
a  calculation  procedure  for  transonic  flow  with  small  amounts  of  trainng- 
edge  separation.  This  method  employs  a  closed  form  solution  of  the  boundary 
layer  equations  in  the  reverse  flow  region.  Limited  comparisons  with  experi¬ 
ment  show  excellent  agreement.  Experience  with  a  similar  procedure  developed 
by  Cousteix  (Ref.  4)  suggests  that  this  approach  can  be  successful  for  air¬ 
foils  with  small  regions  of  separated  flow,  but  must  break  down  at  large 
angles  of  attack  where  wake  modelling  becomes  important. 


(1)  Hicks,  R. ,  Private  communication,  1977. 

(2)  Deiwert,  G.S.,  "Computation  of  Separated  Transonic  Turbulent  Flows", 

AIAA  Paper  75-829,  June  1975. 

(3)  Barnwell,  R.W.,  "A  Potential  Flow/Boundary  Layer  Method  for  Calculating 
Subsonic  and  Transonic  Airfoil  Flow  with  Trailing-Edge  Separation", 

NASA  TM-81850,  June  1981. 

(4)  Cousteix,  J.,  AFOSR-HTTM-Stanford  Conference  on:  Complex  Turbulent  Flows, 
Comparison  of  Computation  and  Experiment,  Volume  I,  Stanford  University, 
SeDtember  1981. 


The  purpose  of  the  present  work  is  the  development  of  an  analysis  method 
for  predicting  the  performance  of  two-dimensional  airfoils  in  transonic  flow 
with  the  main  emphasis  placed  on  its  modelling  of  the  separated  region.  It  is 
important  to  note  that  a  similar  method  for  incompressible  flow,  CLMAX, 
developed  under  Army  Research  Office  supDort  (Contract  DAAG29-76-C-0019,  Ref.  5) 
has  been  very  successful  in  predicting  the  performance  of  two-dimensional  air¬ 
foils  having  large  regions  of  separated  flow.  In  CLMAX,  the  wake  surface  is 
represented  by  a  constant-strength  vortex  sheet  which,  through  iteration, 
assumes  a  force-free  wake  Dosition.  In  the  calculation  method  presented  here, 
TRANMAX,  the  wake  surface  is  modelled  by  a  discontinuity  sheet  with  constant 

-£■$  along  it  and  with  a  jump  in  <J>  (tangential  velocity)  across  it.  This 

d$ 

approach  is  analogous  to  the  constant  strength  vorticity  sheet  model,  CLMAX. 

The  wake  model  in  the  present  code  is  closed  and  the  shape  remains  fixed  through 
one  complete  inviscid  flow  iteration  cycle. 

Jameson's  code,  FL06  (Ref.  6),  with  substantial  change  to  include  the 
wake  region,  is  used  for  the  inviscid  flow  calculation.  The  Cohen-Reshotko 
(Ref.  7)  laminar  boundary  layer  method  and  Green's  (Ref.  8)  turbulent  lag- 
entrainment  boundary  layer  method  are  employed  for  the  viscous  flow  calcula¬ 
tions.  The  details  of  the  procedure  and  its  performance  are  described  in 
subseguent  sections. 


(5)  Dvorak,  F.A.,  Maskew,  B.  and  Rao,  B.M. ,  "An  Investigation  of  Separation 
Models  for  the  Prediction  of  Co  ",  Final  Technical  Report,  Contract 
DAAG29-76-C-0019,  Prepared  for  m  the  U.S.  Army  Research  Office, 

Research  Triangle  Park,  N.C.,  April  1979. 

(6)  Jameson,  A.,  "Numerical  Computation  of  Transonic  Flows  with  Shock  Waves", 
International  Uni  on  of  Theoreti cal  and  Applied  Mechanics ,  Springer-Verlag 
New  York,  Inc.,  September  1975,  pp.  384-414. 

(7)  Brune,  G.W.  and  Manke,  J.W.,  "An  Improved  Version  of  the  NASA  Lockheed 
Multielement  Airfoil  Analysis  Computer  Program",  NASA  CR-15323,  March  1978, 
pp.  69-87. 

(8)  Green,  J.E. ,  Weeks,  D.J.  and  Brooman,  J.W.F.,  "Prediction  of  Turbulent 
Boundary  Layers  and  Wakes  In  Compressible  Flow  by  a  Lag-Entrainment  Method", 
Royal  Aircraft  Establishment  TR-72231  ,  December  1972. 


2.0  DESCRIPTION  OF  THE  ANALYSIS  METHOD 

2.1  General  Description 

The  present  analysis  method  is  built  upon  Jameson's  two-dimensional 
potential  flow  calculation  method  (Ref.  6),  FL06,  Cohen/Reshotko's  laminar 
boundary  layer  method  (Ref.  7)  and  Green's  lag-entrainment  turbulent  boundary 
layer  method  (Ref.  8).  These  methods  were  verified  thoroughly  by  their 
originators  and  can  be  used  for  transonic  flow  calculations  with  confidence. 

The  solution  procedure  is  shown  in  Figure  1.  The  initial  potential  flow 
solution  is  obtained  either  with  or  without  the  wake  prescribed,  depending 
on  the  initial  seoaration  point.  The  potential  flow  method  has  been  sub¬ 
stantially  modified  to  model  the  separation  region  with  the  shape  of  the 
wake  generated  by  a  procedure  to  be  described  in  a  later  section.  Having 
obtained  the  inviscid  pressure  distribution,  the  boundary  layer  development 
is  predicted  using  the  methods  mentioned  earlier.  This  completes  one  viscid/ 
inviscid  calculation  cycle. 

In  subseguent  iterations,  the  induced  normal  velocity  on  the  airfoil 
surface  due  to  boundary  layer  displacement  effect  is  taken  into  account.  A 
new  wake  is  generated  at  the  start  of  every  iteration. 

This  procedure  is  repeated  until  the  solution  converges;  i.e.,  the 
separation  points  between  two  successive  iterations  remain  unchanged. 

Details  of  individual  elements  are  fully  described  in  the  following 
sections. 


2.2  Potential  Flow  Calculation  Method 
2.2.1  Basic  Eguations 

From  the  eguation  of  continuity  and  the  momentum  eguation,  we  have 
<a2-u2)f  -uv(f +  f)  +  (a*-v2)f  =  0  (1) 


Assuming  irrotational  flow, 


a  velocity  potential,  <f>,  can  be  defined;  i.e.. 


u 


3x  ’ 


(2) 
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Figure  1.  Calculation  Procedure. 


Substitute  these  into  Eqn.  (1)  to  obtain: 

(a2  -  u2)$xx  -  2uv<j>xy  +  (a2  -  v2)4>yy  =  0  (3) 

This  equation  can  be  solved  for  <j>  with  the  use  O'f  the  energy  equation 

a!  ■  (^H1)  q-!  '4> 


where 

q  =  y'u2  +  v2 


The  Neumann  boundary  condition  is  prescribed  along  the  surface:  =  0,  is 

oP 

set  initially,  and  it  takes  a  new  value  which  reflects  the  viscous  effect  after 
each  iteration  cycle 


=  0  (First  Iteration) 


^ -^-(pufi*)  =  f  (Subsequent  Iterations) 


(5) 


where  <5*  is  the  boundary  layer  displacement  thickness  and  s  is  measured  alonq 
the  airfoil  surface. 

Because  of  the  nature  of  Eqn.  (3),  hyperbolic  if  the  local  Mach  number, 

M  =  q/a  >  1  and  elliptic  if  M  <  1,  it  is  essential  to  transform  the  infinite 
flow  field  onto  a  finite  domain.  This  can  be  achieved  by  mapping  the  exterior 
of  the  airfoil  in  the  z-plane  conformally  onto  the  interior  of  a  unit  circle 
in  the  cr-olane  as  shown  in  the  following  figure  taken  from  Reference  9. 


(9)  Roqers,  E.O.,  "Numerical  Solution  of  Subcritical  Flow  Past  Airfoils".  NSRDC 
Reoort  4112,  May  1973. 


Figure  2.  Grid  System. 


This  transformation  is  particularly  useful  because  an  evenly  distributed 
qrid  system  in  the  circle  plane  gives  denser  finite  mesh  near  the  body  at 
the  leading  and  tracing  edges  in  the  physical  plane  where  it  is  needed  most. 


2.2.2  Mapped  Coordinate  System 
Consider  the  speed  in  the  z-plane. 


where 


Thus  we  have 


=  -M./R 
ax  rae/t5’ 


M  =  li/B 
ay  ar/D 


(7) 
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The  complex  potential  about  a  unit  circle  in  a  uniform  stream  is  given 


by: 


w  =  UM(z  +  j)  +  iEJlnz 


(8) 


Here  the  velocity  potential  at  infinity  becomes  unbounded  and  is  of  the  order 
R  Since  the  present  transformation  requires  the  direct  inversion  of  this 
external  flow,  the  singular  behavior  at  infinity  inevitably  occurs  at  the 
origin  of  the  a-plane. 

In  order  to  remove  this  singularity  at  the  origin,  0(l/r),  and  the  dis¬ 
continuity  at  9  =  2-tt  due  to  the  circulation,  a  translated  potential,  G,  is 
i  ntroduced. 


G  = 


+  E(e  +  o) 


(9) 


where  2ttE  is  the  circulation,  and  a  is  the  angle  of  attack. 

Substituting  Eq ns.  (7)  and  (9)  into  Eqn.  (2)  to  obtain 

(a2  -  u2)  G00  -  2uvrGr0  +  (a2  -  v2)  r  ~  (rGr) 

-  2uv(Gq  -  E)  +  (u2  -  v2)Gf  +  (u2  +  v2)(^  B0  +  vBr)  =  0  (10) 

where 

r(GQ  -  E)  -  sin  (0  +  a)  r2G  -  cos  (0  +  a) 

..  ^  ^  ... _ \i  =  -  —  ■  ■  '  ■ 


The  Neumann  boundary  condition  reduces  to 


Gf  =  cos  (8  +  a) 


f 

B 


at  r  =  1 


(ID 
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while  the  far-field  boundary  condition  becomes 


G 


6  +  a  -  tan 


at  r  =  0 


(12) 


Here,  the  circulation  constant,  E,  is  determined  by  the  Kutta  condition 
derived  from  the  upper  and  lower  surface  separation  pressures  and  is  discussed 
in  a  later  section. 


2.2.3  Inviscid  Model  of  Separated  Region 

As  mentioned  earlier,  the  modelling  of  the  wake  is  one  of  the  most 
important  elements  in  the  method.  One  may  solve  the  parabolized  Navier  Stokes 
equations  for  the  special  class  of  flow  instead  of  modelling  the  wake.  However, 
this,  in  principle,  cannot  handle  the  region  of  large  reverse  flow.  On 
the  other  hand,  solving  the  full  time-dependent  Navier  Stokes  equations  is 
not  practical  at  this  time  because  of  the  enormous  computing  time  required. 

The  only  logical  approach  to  this  oroblem  is  to  simulate  the  separated  flow 
with  an  inviscid  wake  model.  A  good  wake  model,  in  fact,  can  represent  the 
actual  flow  very  well  as  was  demonstrated  by  an  earlier  study  (CLMAX). 

In  the  present  model,  the  separated  region  is  confined  by  two  dividing 
streamlines,  one  from  the  upper  surface  separation  point  and  the  other  from 
the  trail inq-edqe  point,  as  shown  in  Figure  3. 


Figure  3.  Wake  Model. 

As  shown  in  the  figure,  along  the  dividing  streamlines  the  gradient  in  velocity 
potential,  ||,  is  constant,  while  across  the  streamlines  |£  is  zero. 
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The  wake  shaoe  must  be  updated  according  to  the  new  separation  point. 

If  the  separation  point  remains  unchanged,  then  the  solution  is  considered 
to  be  converged. 

Construction  and  Transformation  of  the  Wake 

Once  the  separation  point,  S,  on  the  upper  surface  is  found,  the  wake 
is  generated  in  the  following  manner  (see  Fiqure  3).  M  is  the  mid-point  on 
the  chord,  ST,  and  the  point,  J,  i_s  fixed  at  twice  the  wake  width  (Wo)  away 
from  the  Doint,  M.  The  slope  of  MN  is  chosen  to  be  .45  times  the  trailing- 
edqe  sloDe.  Two  parabolas  can  then  be  constructed  through  SN  and  TN,  whose 
initial  sloDes  are  equal  to  the  body  slopes  at  their  respective  points. 

This  wake  in  the  physical  Diane  can  now  be  transformed  onto  the  circle 
Diane.  First,  the  grid  network  on  the  physical  plane,  which  corresponds  to 
the  uniform  grid  in  the  circle  plane,  is  generated  by  usinq  the  mapping 
function.  With  the  aid  of  simple  interpolation,  the  wake  in  the  physical 
plane  can  be  readily  transformed  onto  the  circle  plane.  A  typical  resulting 
separated  region  in  the  circle  plane  is  shown  in  Fiqure  4. 


Fiqure  4.  Separated  Region  in  the  Computational  Plane. 

The  streamlines  which  divide  the  flow  into  attached  and  separated  flow 
regions  represent  lines  of  constant  ||  .  As  in  the  wake  model  used  in  Program 

CLMAX,  these  streamlines  represent  thin  shear  layers  across  which  there  exists  a 
a  lump  in  tanqential  velocity;  and,  consequently,  a  jump  in  velocity  potential. 
Due  to  this  jump,  a  special  treatment  in  the  finite-difference  scheme  is 
necessary  along  this  boundary.  In  the  present  scheme,  the  value  of  4>  is  adjusted 
up  or  down  by  the  amount  of  jump,  A<J>,  where 

A<f)i  =  H  Asi 
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and  as  is  measured  from  point  S  or  T,  depending  on  the  particular  side,  to 
make  the  function  <j>  continuous  across  the  dividing  streamline.  With  this 
addition  and  subtraction,  the  same  difference  formulae  can  be  used  throughout 
the  field.  The  procedure  is  similar  to  the  so-called  "shock-fitting"  technique. 

The  Kutta  condition,  which  is  normally  applied  at  the  trailing  edge,  is 
aoplied  at  points  S  and  T;  i.e.. 


3s 


S 


(13) 


Special  Treatment  in  Finite-Difference  Formulae  Along  the  Dividing  Streamline 


=  Constant 


Figure  5.  Finite-Difference  Grid. 

The  main  idea  of  present  finite-difference  formulae  near  the  dividing 
streamline  is  to  make  the  potential,  <j>,  continuous  across  this  streamline. 
As  previously  mentioned,  this  can  be  achieved  by  adding  or  subtracting  the 
amount  of  jump,  A<f>,  to  or  from  <p  a*,  a  particular  point.  The  correction  Is 
applied  to  all  points  lying  on  the  side  opposite  the  point  of  finite-dif¬ 
ference  approximation. 
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As  for  examples,  the  derivatives  about  the  point,  0,  are  shown  below. 


a<j>/ax  = 

(4>r  -  <t>j_ ) / ( 2  •  Ax) 

(14(a)) 

84>/aY  = 

(*p  -  <*„  -  A*n))/(2  ■  AY) 

(14(b)) 

where  Ae„  ■  ff I  •  As„ 

at  separation 

As^  =  arc  lenqth  along  the  dividing  streamline 


324>/3X2  =  (<^  -  2$q  +  <i>^)/Ax2 

(14(c)) 

32<t/3Y2  =  | p  -  2^q  +  (<J>^  -  A4>n ] 

l)/AY2 

(14(d)) 

32<j)/axaY  =  +  i'K  -  Cd>q  + 

A*Q))/4(AXA  Y) 

(14(e)) 

2.2. 4  Method  of  Solution 

The  transformed  Eqn.  (10)  with  boundary  conditions,  Eqns.  (11) 
and  (12),  is  solved  by  a  finite-difference  scheme.  Upwind  differencing  is 
used  when  the  local  flow  is  supersonic  ,  while  central  difference  formulae 
are  used  at  subsonic  Doints.  The  resulting  set  of  difference  equations  is 
solved  iteratively.  Details  of  the  method  are  referred  to  in  the  original 
paper  (Ref.  6)  and  only  the  major  differences  are  presented  here. 

The  Kutta  condition  is  imposed  by  matching  the  pressures  at  the  upper 
surface  separation  point  and  the  trailing  edge;  i.e.. 


as 


sep 


-  M. 

'  3s 


trailing  edge 


(13) 
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The  amount  of  circulation,  E,  in  Eqn.  (9)  is  continuously  updated  to  satisfy 
Eqn.  (13).  The  new  estimate  of  E,  determined  from  the  Kutta  condition,  is 
then  used  for  the  next  iteration.  This  cycle  is  repeated  until  the  correction 
to  E  satisfies  the  converqence  criteria. 

It  was  found,  during  the  course  of  this  work,  that  the  method  may  encounter 
some  difficulty  if  the  solution  is  sought  directly  for  free  stream  Mach  numbers 
greater  than  0.35.  This  occurs  because  the  actual  flow  is  becomming  less 
stable  and  will  respond  rapidly  to  even  small  disturbances.  When  the  separation 
region  is  large  (greater  than  10%  chord),  the  initial  inviscid  flow  field  is 
not  representative  of  the  flow  field  with  separation;  consequently,  a  poor 
initial  solution  can  lead  to  divergent  behavior  at  higher  Mach  numbers.  Proper 
converqence  and  stability  can  be  achieved  by  starting  the  solution  from  a  low 
Mach  number  and  raising  the  Mach  number  incrementally  to  the  desired  value. 

After  the  translated  potential,  G,  is  obtained,  the  tanqential  velocity 
component,  u,  on  the  surface  can  be  obtained  readily  throuqh  simple  dif¬ 
ferencing.  The  pressure  coefficient,  C  ,  along  the  surface  is  given  by: 


2.3  Calculation  of  Boundary  Laver  Development 

As  in  the  potential  flow  calculation  procedure,  new  methods  for  both 
laminar  and  turbulent  boundary  layers  were  adopted.  Curie's  method  for 
laminar  flow  and  the  Nash-Hicks  method  for  turbulent  flow  are  replaced  by  the 
Cohen-Res hotko  method  and  Green's  method,  respectively.  However,  the  same 
transition  criterion,  the  Granville  procedure,  is  still  intact.  The  possi¬ 
bility  of  the  reattachment  as  a  turbulent  boundary  layer  after  laminar  boundary 
layer  separation  is  also  examined  on  the  basis  of  the  Reynolds  number  based 
on  the  momentum  thickness  at  the  point  of  separation  (Ref.  10). 


2.3.1  Cohen  and  Reshotko  Method 

This  integral  method  was  developed  for  the  steady  two-dimensional 
laminar  boundary  layer  with  the  assumption  that  the  surface  temperature  is 
uniform. 


(10)  Dvorak,  F. A.  and  Woodward,  F.A.,  "A  Viscous/Potential  Flow  Interaction 
Analysis  Method  for  Multi-Element  Infinite  Swept  Winqs;  Volume  I",  NASA 
CR-2476,  November  1974. 
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First,  consider  the  Stewartson  transformation. 


3  it)  pu 

ay  “  p0  ‘ 


3^  _  -pv 
ax  "  p0 


x 


X 


o 


-£-  dY 


u 


ay 


(16) 


where  ^  and  a  represent  the  stream  function  and  the  speed  of  sound,  respectively. 
Capital  letters  denote  quantities  of  the  equivalent  incompressible  problem. 

With  this  transformation,  the  equivalent  incompressible  problem  can  be  formu¬ 
lated.  Now,  define  a  parameter, 


n 


e’tr  dUe 


(17) 


where  subscript,  tr,  denotes  the  transformed  coordinate.  After  some 
manipulation,  one  can  obtain 


n  -  -  ClMe 


-C2^e 
dx 


(18) 


where  and  C2  are  constants. 

This  can  be  readily  solved  for  n  by  simple  integration  formulae  such  as 
the  trapezoidal  rule.  Then  the  momentum  thickness  and  the  shape  factor  are 
obtained  from  the  explicit  expressions  of  n.  (See  Ref.  7  for  details.) 


The  separation  of  the  boundary  layer  is  detected  by  examining  the 
Pohl  hausen  parameter,  *  assuming  that  the  H-A  table  for  the  in¬ 

compressible  flow  is  still  valid. 

After  transition,  either  through  natural  transition  or  through  a  sepa¬ 
ration/reattachment  process.  Green's  method  takes  over  the  calculation  of  down¬ 
stream  turbulent  boundary  layer  development. 


2.3.2  Green's  Method 

This  is  a  "lag-entrainment"  integral  method  involving  three  equations 
momentum- integral ,  entrainment  and  a  rate  equation  for  the  entrainment  coeffi¬ 
cient  (Ref.  8).  The  method  is  a  combination  of  Head's  original  entraiment 
momentum-integral  method  and  the  turbulent  model  proposed  by  Bradshaw  et  al. 
in  which  the  algebraic  relation  for  the  entrainment  coefficient  of  Head's 
method  is  replaced  by  a  rate  equation  derived  from  the  turbulent  kinetic  energy 
equation. 

The  resulting  equations  are: 


CE> 

(19) 

(20) 

dCr 

W  =  H3  (9,  H’ 

(21) 

where  » 

•’/ 

o 

00 

*  - 1  /  f  <’  -  f  >  oy 

9  J  pe  e 
o 

0 
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This  system  of  ordinary  differential  equations  is  solved  by  the  Runge-Kutta 
method.  The  separation  point  can  be  located  by  monitoring  both  the  friction 
coefficient,  C^,  and  the  shape  factor,  H. 


3.0  RESULTS  AND  DISCUSSION 


The  main  purpose  of  this  analysis  method,  TRANMX,  is  to  predict  the 
performance  of  a  given  airfoil  for  a  wide  range  of  angles  of  attack,  and  to 
determine  the  maximum  lift  and  its  corresDonding  angle.  As  will  be  presented 
in  this  section,  TRANMX  performs  very  well  for  all  the  cases  examined.  The 
oressure  distributions  and  the  C.  -  a  curves  are  in  excellent  agreement  with 
the  experimental  data. 

The  method  has  been  tested  against  four  distinctly  different  airfoils; 
i.e.,  GA(W)-1,  NACA  0012,  Model  A1  and  a  Huqhes  Helicopter  Comoany  airfoil 
with  tab.  The  free  stream  Mach  number  varies  from  .15  for  the  GA(W)-1  to 
.6  for  the  A1  airfoil.  For  each  flow  condition,  the  angle  of  attack  is 
gradually  increased  until  the  airfoil  has  stalled.  Airfoils  at  angles  of 
attack  beyond  the  static  stall  angle  have  been  analysed. 

Figure  6(a),  (b)  and  (c)  shows  the  pressure  distributions  about  the 
airfoil,  GA(W)-1,  for  M^  =  0.15  and  Re  =  6  x  106.  It  should  be  noted  that 
this  is  the  same  test  case  with  which  the  results  of  CLMAX  were  in  good  agree¬ 
ment.  The  good  correlation  shown  in  these  figures  confirms  that  these  two 
analytical  methods  are  compatible.  The  computed  is  a  little  higher  than  the 
data  as  shown  in  Figure  7  due  to  the  slightly  higher  pressure  distribution 
along  the  lower  surface.  However,  this  discrepancy  of  less  than  3%  is  well 
within  the  experimental  accuracy. 

The  comparisons  shown  in  Figures  8  and  9  are  for  the  NACA  0012  symmetric 
airfoil  at  M^  =  0.5,  Re  =  3  x  10  ,  Here  again  the  results  are  quite  good 
although  the  quoted  angles  of  attack  appear  to  be  high.  It  is  interesting 
to  note  that  the  pure  inviscid  solution  for  a  =  9.86,  which  is  plotted  in 
Figure  8(c),  fails  to  predict  the  real  pressure  distribution  by  a  surprising 
margin.  A  separated  region  of  about  20%  of  the  chord  makes  the  pure  poten¬ 
tial  flow  solution  meaningless.  This  indicates  that  the  viscid/inviscid 
interation  changes  the  real  flow  field  considerably  and  the  viscous  effects 
cannot  be  ignored. 

The  results  of  two  different  calculation  methods  are  plotted  in  Figure  10 
along  with  experimental  data  for  a  Hughes  Helicopter  Company  airfoil.  The 
free  stream  Mach  number  is  .46  and  the  Reynolds  number  is  3.8  million.  CLMAX 
does  not  perform  very  well,  particularly  because  it  is  not  developed  to  handle 
the  highly  compressible  flow.  This  also  shows  the  limitation  of  the  usage 
of  linearized  equations  in  compressible  flow  for  this  Mach  number.  TRANMAX 
performs  well  except  at  the  highest  angle  of  attack  (Figure  10(c)).  In  this 
case,  the  actual  flow  reaches  a  local  Mach  number  of  1.5  and  the  boundary 
layer  is  believed  to  separate  at  the  foot  of  the  shock.  This  boundary  layer 
(separation)  effectively  removes  the  distinct  shock  but  fails  to  reattach. 

This  separated  flow  exhibits  quite  different  characteristics,  e.g.,  varying 
pressure  along  the  body,  and  cannot  be  treated  as  well  by  the  present  method. 
The  Co  -  a  curve  shown  in  Figure  11  is,  therefore,  rather  fortuitous  where 
a  >  10°. 
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(a) 

Figure  6.  Pressure  Distribution  along  GA(W)-1  Airfoil. 
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(b) 

Flqure  6.  Continued. 
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Fiqure  6.  Concluded. 
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qure  7.  -  a  Curve  for  GA(W)-1  Airfoil. 
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Flqure  8.  Concluded. 
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(a) 

Figure  10.  Pressure  Distribution  along  a  Hughes  Helicopters 
Company  Airfoil. 
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Another  helicopter  section,  the  A1  airfoil,  results  are  compa1'  ‘  in 
Fiqures  12,  13  and  14.  The  data  was  collected  by  Hicks  and  McCroskey  (Ref.  11) 
in  the  Ames  2'  x  2'  Transonic  Wind  Tunnel  with  the  model  airfoil  chord  of  0.5 
feet.  Among  a  wide  ranqe  of  test  conditions,  only  the  Mach  numbers  of  0.5 
and  0.6  are  chosen  owing  to  the  present  interest.  The  Cj,  -  a  curve  shown 
in  Figure  12  can  be  divided  into  three  reqions.  The  flow  in  the  first  region, 
a  <  8°,  can  be  characterized  by  the  tr ailing-edge  separation.  The  lift,  C£, 
varies  nearly  linearly  with  respect  to  a  in  this  region.  The  next  region 
represents  the  flow  for  anqles  of  attack  from  8°  to  aCLMAX*  boundary 
layer  in  this  region  separates  at  the  foot  of  the  shock  and  thus  blunts  the 
shock  front.  In  this  region  a  combination  of  shock  separation  and  trai ling- 
edge  separation  are  present.  Cjj.  varies  gradually  in  this  region  until  the 
catastrophic  leading-edge  separation  brinqs  the  lift  down  sharply  in  the 
third  reqion. 

The  results  are  in  good  agreement  with  experiment  in  region  I.  The 
calculated  is  much  higher  than  the  experimental  data  in  the  region  II, 
although  the  general  trend  of  C&  -  ot  curve  is  in  close  agreement  with  the 
experiment,  the  lower  experimental  C&  is  due  to  the  shock  boundary  layer 
separation  which  the  present  method  currently  cannot  predict.  The  pressure 
distributions  shown  in  Figures  13  and  14  correlate  very  well  with  the 
experimental  data.  It  is  interesting  to  note,  in  Figure  14,  the  unsteady 
behavior  of  the  pressure  in  the  leading-edge  region.  This  unsteadiness  may 
indicate  the  temporary  presence  of  shock  boundary  layer  separation  bubble  which 
was  mentioned  earlier. 

As  we  have  seen  so  far  in  this  section,  the  new  calculation  method  performs 
very  well  for  all  the  cases  tested.  The  minor  discrepancies  which  exist  in 
some  of  the  conroarisons  are  attributed  to  the  calculation  method's  own 
limitations  rather  than  inconsistent  performance.  Misleading  high  lift,  which 
is  likely  to  happen  when  the  method  does  not  detect  the  shock  boundary  layer 
separation,  can  be  easily  noted  by  a  careful  analysis  of  the  result.  One 
obvious  clue  may  be  the  unrealistically  high  local  speed,  say,  M  >  1.4,  since 
the  boundary  layer  through  this  shock  is  not  likely  to  remain  attached. 

The  main  difficulty  in  verifying  the  calculation  method  has  been  the  lack 
of  quality  transonic  flow  data.  Most  of  the  available  data  are  limited  to  the 
case  of  no  or  small  seDarated  flow  regions.  In  some  cases,  e.g.,  the  A1  airfoil, 
the  angle-of-attack  correction  due  to  tunnel  blockage  effect  is  so  large  (for 
M  =  0.6,  Aa  ~  a)  that  the  credibility  of  the  data  is  in  serious  doubt.  More  . 
reliable  high  Mach  number  flow  data  with  a  substantial  reaion  of  reversed 
flow  are  desired  for  the  future  improvement  of  the  calculation  method. 


(11)  Hicks,  R.M.  and  McCroskey,  W.J.,  "An  Experimental  Evaluation  of  a  Heli¬ 
copter  Rotor  Section  Designed  by  Numerical  Optimization" ,  NASA  TM-78622, 
March  1980. 
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Fiqure  13.  Pressure  Distribution  along  an  A1  Airfoil 
for  =  0.596;  a  =  6°. 
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4.0  CONCLUSIONS 


An  analysis  method  for  transonic  trai ling-edge  type  separated  flow  has  been 
successfully  demonstrated  in  comparison  with  experimental  data.  Predicted 
Ca  -  a  behavior  is  in  good  agreement  with  experiment  for  several  different  air¬ 
foil  types  at  Mach  numbers  approaching  0.6.  At  higher  Mach  numbers  inspection 
of  the  local  Mach  number  distribution  usually  suggests  separation  as  a  result  of 
strong  shocks.  While  the  program  is  not  currently  capable  of  modelling  this 
phenomenon,  it  can  be  used  as  a  predictive  tool  to  indicate  when  shock/boundary 
layer  separation  can  be  expected. 

The  main  advantage  of  this  method  over  other  existing  methods  is  that  the 
present  scheme  with  its  simple  wake  model  allows  computations  for  massive  sepa¬ 
ration.  It  also  can  be  extended  to  the  three-dimensional  case  where,  as  pre¬ 
viously  cited,  alternate  methods  are  limited  to  two-dimensional  flow. 
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